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I. INTRODUCTION 

There are several reasons why polarized M0ller scattering has attracted so much interest from both experi- 
mental and theoretical communities. The scattering of two identical polarized fermions has been used in many 
laboratories for the high-precision determination of the electron-beam polarization, including SLC in [lj], E-143 
and E-154 0] experiments at SLAC, and several experiments at JLab Q and MIT-Bates @. A M0ller po- 
larimeter may also be of use in future experiments planned at ILC @. An experiment E-158 at SLAC Q, which 
studied M0ller scattering of 45- to 48-GcV polarized electrons on unpolarized electrons of a hydrogen target, 
allowed one of the important parameters in the Standard Model - sine of the Weinberg angle - to be determined 
with an unprecedented accuracy. The new-generation experiment at 11 GeV soon to be started at JLab will 
offer a new level of sensitivity and measure the parity-violating asymmetry in the scattering of longitudinally 
polarized electrons off unpolarized electrons to a precision of 0.73 ppb. That would allow a measurement of the 
weak charge of the electron to a fractional accuracy of 2.3% and a determination of the weak mixing angle with 
an uncertainty of ± 0.00026 (stat.) ± 0.00013 (syst.) 0. 

Since e~e _ — > e~e~("/) scattering is a very clean process with well-known initial energy and extremely 
suppressed backgrounds, any inconsistency with the Standard Model will signal new physics. M0llcr scattering 
experiments can provide indirect access to physics at multi-TeV scales and play important complementary role 
to the LHC research program |{|. 

Obviously, before we can extract reliable information from experimental data, it is necessary to take into 
account higher order effects. These are the processes which are more complicated than the process being 
studied, but which are indistinguishable from it experimentally. This procedure of the inclusion of radiative 
corrections is an indispensable part of any modern experiment, but will be of the most paramount importance 
for the ultra- precise measurement of the weak mixing angle via 11 GeV M0ller scattering planned at JLab. 

One of the earliest works on Electroweak Radiative Corrections (EWC) to the observables of M0ller scattering 
was done by Czarnecki and Marciano According to their calculations for the asymmetry in the E-158 

kincmatical region, the EWC reduce the tree level prediction by 40±3%. A close value of the asymmetry was 
obtained in the study of Denner and Pozzorini (Tl| , where radiative corrections in polarized M0ller scattering 
were examined at arbitrary, including high, energies. It is worth noticing that these two studies used different 
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renormalization schemes: [13] employed the modified minimal-subtraction scheme and 11 1 used the on-shell 



rcnormalization. However, the authors of [13] and [HI did not take into account all radiative contributions 
consistently. For example, they completely disregarded the Hard Photon Bremsstrahlung (HPB) contribution 
assuming that it was small. The first computation of the HPB contributions with realistic experimental cuts was 
done by Petriello [l2[ who studied a total set of the lowest-order EWC under kinematical conditions of E-158 
experiment. He clearly showed that experimental cuts play a significant role. Later works [l3|, [lH and [l5[ 
dedicated to the EWC for the E-158 experiment which employed a covariant method for extracting of infrared 
divergences [l6j showed a good agreement with [Io| for the bulk of weak-interaction contributions and proved 
that the effect of HPB was below the experimental error of E-158. 

Electromagnetic radiative corrections, which dominate weak interaction effects at the energies of E-158, and 



their effect on the measurement of beam polarization in a M0ller polarimeter were assessed in 17j . A different 
calculation of these electromagnetic corrections was undertaken in [l8| j where attention was primarily given to 
the effect of experimental cuts on the inelasticity and, accordingly, proving the need for taking into account 
hard-photon bremsstrahlung in experiments of this kind. A Monte Carlo generator that makes it possible to 
simulate radiative events in a M0ller process at moderate energies with the aim of determining beam polarization 
was developed in [l9| . 

Finally, an attempt to calculate the hard bremsstrahlung contribution for M0ller process at high energies was 
made in paper [23] ■ A feature specific to this calculation, as well as to the calculations in [ll[ and [l2[ , was the 
use of a parameter u) that separates regions of soft and hard photons. As a result, the Monte Carlo integration 
technique permitted the development of a code which made it possible to take into account radiative corrections 
in the case of arbitrary cuts imposed on the detection of electrons. 

Although obviously significant theoretical effort has been already dedicated to this task, we can see at least 
three major reasons for addressing EWC to M0ller scattering yet again. First, since the new experiment at 
Jlab will provide more precise information than, for example, E-158, the theoretical predictions for this ultra- 
precise measurement must be made with a full treatment of one-loop radiative corrections and at least leading 
two-loop corrections. In this work, we calculated the full set of the one-loop corrections both numerically with 
no simplifications using FeynArts and FormCalc [2l| as the base languages and analytically in a compact form 
useful for fast estimations. This way, we can control an accuracy of asymptotic approximations at the one-loop 
level very well now, so we hope that the same idea will help us to control our accuracy as we calculate leading 
corrections at two-loop level in the nearest future. The next major goal of the presented paper is to get EWC 
in a form which is analytically free from nonphysical parameters. And, finally, we believe that the complexity 
of the problem demands the tuned comparison with different calculation schemes. This will be our next task. 

The rest of the paper is organized as follows. The basic notations as well as the lowest-order (Born) contri- 
bution to Moller scattering are presented in Section II. The contribution from a full set of additional virtual 
particles (such as self-energies, vertices and boxes) is described in Section III. The real photon emission and its 
separation on hard and soft parts as well as cancellation of the infrared divergences together with the unphysical 
parameters can be found in Section IV. The numerical analysis is presented in Section V and Appendix B. 

II. BORN CROSS SECTION: BASIC NOTATION AND CONVENTIONS 

In the Standard Model, the Born cross section for the M0ller scattering 

e"(fei,£) + e~{jp 1 ,ifj e~(k 2 ) + e~(p 2 ) (1) 

can be represented in the form 
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cr° = V [>&(u 2 D it D jt + t 2 D iu D lu ) + A^'s 2 (D lt + D lu )(D ]t + D ju )], (2) 

s L — ' 



with a = da jd cos 9 where 9 is the scattering angle of the detected electron with the 4-momentum k 2 in the 
center of mass system of the initial electrons and the 4- momenta of initial (fci and p\) and final (k 2 and p 2 ) 
electrons (see Fig. [T]) generate a standard set of Mandelstam variables, 

s = (k 1 +p 1 ) 2 , t=(fci-fc 2 ) 2 , u={k 2 - Pl f, (3) 

while the electron polarization vectors £ and r\ have the form [22j] : 

k x 2m 2m p x 

£~ Pi, k i • ( 4 ) 

m s s m 

It should also be noted that the electron mass m is disregarded wherever possible, in particular if m 2 -c s, —t, —it. 



(1) 



(2) 



FIG. 1: Diagrams describing nonradiative M0ller scattering in the (1) t- and (2) u-channels. 



A useful structure in the present study is 

£>« = -J-j (i=7,Z; l = t,u), (5) 

which depends on the Z-boson mass mz and on the photon mass m 7 . The photon mass is set to zero everywhere 
with the exception of specially indicated cases where the photon mass is taken to be an infinitesimal parameter 
that regularizes an infrared divergence. Another set of useful functions is 

\ i.k \ i,k \ i.k , \ i,k \ i,k / o\ 

a± —Mb Mt ± a 2b Mt i y°) 

which are combinations of coupling constants and Pb(t) ; where Pb(t) are the degrees of polarizations of electrons 
with 4-momentum k\ (pi). Specifically, they are given by 

M B (T) ~ "V -PB(T)A A i 2 B(T) ~ A A ~ PB(T)Ay , 

X l y = v i v j + aV, X\ j = v i a j + aV, (7) 

where 

^ = 1, a? = 0, v z = {II + 2s w )/(2s w c w ), a z = I^/(2 Sw c w ). (8) 

It should be recalled that J| = —1/2 and sw (cw) are the sine (cosine) of the Weinberg angle expressed in 
terms of the Z- and W-boson masses according to the rules of the Standard Model: 

c w = m w /m z , % = \Jl- Cyy- ( 9 ) 

The electron polarization degrees Pb(t) m the cross sections are labeled here as follows: the subscripts L 
and R on the cross sections correspond to Pb(t) — — 1 an d Pb(T) — +1 5 where the first subscript indicates the 
degree of polarization for the 4-momentum k\, while the second one indicates the degree of polarization for 
the 4-momentum p\. By combining the degrees of electron beam polarizations, we can obtain four measurable 
cross sections, but, by virtue of the rotational invariance, the two of them will be identical: olr = <jrl- From 
the three cross sections we can construct three independent asymmetries [23| (which are very close at large 
scattering angles), and the two of them (mainly A{) are main subject of our investigation: 

^ _ Q"LL + OLR - CTRL - VRR _ PLL - &RR 

&LL + &LR + &RL + &RR O LL + 2(7 LR + O RR ' 

A 2 = a ^- a ™ . (n) 

OLL + <JRR 

All of the asymmetries are proportional to the combination 1 — 4s w (by virtue of the proportionality of the 
cross-section difference ctll — vrr) and are therefore highly sensitive to small changes in sw- This is precisely 
the reason why the asymmetry A±, which, at moderately low energies, is given by 



(12) 
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FIG. 2: One-loop t-channel diagrams for the M0ller process. The circles represent the contributions of self-energies and 
vertex functions. The u-channel diagrams are obtained via the interchange fa Hp2. 



was used as the observable in the E-158 and will be measured in the future experiment at JLab with 11 GeV 
electrons. 

In addition to a fact that the asymmetry A\ is a parity-violating observable, this asymmetry has yet another 
remarkable property based on its structure: it only involves combinations o~ll + &lr and grl + vrr, which 
can be interpreted as the cross sections for the scattering of electrons having the polarizations of ps = —1 and 
Pb = +1 on unpolarized electrons. Frequently, A\ is said to be a single-polarization asymmetry, in contrast to 
the asymmetry constructed as 

A LR = ^7 LL , (13) 

&LR + &LL 

which consists of the cross sections for the scattering of electrons having different polarizations. The latter asym- 
metry, which conserves parity, is the most important observable in determining the electron-beam polarization 
with the aid of a M0llcr polarimeter. 



III. EWC: CONTRIBUTION OF ADDITIONAL VIRTUAL PARTICLES 



A. Boson Self-Energies 

The contribution of virtual particles (^-contribution) to the observables of Moller scattering is described by 
three classes of diagrams (see Fig. boson self-energies (BSE), vertex functions, and two-boson exchange 
(boxes). In the on-shell renormalization schemes which we use here there are no contributions from the electron 
self-energies. 

The corresponding cross section is given by 

a v = a BSE + ( 7 VCT + a Box , (14) 

where the first term corresponds to the contributions of the boson self-energies, the second term represents the 
vertex diagrams, and the third term stands for the diagrams of two-boson exchange. 

The contributions of the photon and Z-boson self-energies are shown symbolically in Fig. [UJl). They do not 
involve infrared divergences and have the following form 

- BSE = ^£A/f + (^), (15) 
l=i 

where are expressed in terms of the propagators and functions M ev ,od 

Mli kl = 2{s 2 +u 2 )\ 1 %\ 1 k l + 2{s 2 -u 2 )\ 2 %\ 2 k l, 

M*f = -2s 2 (X 1 %X 1 k J+X 2 %X 2 k 1 i), (16) 
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are defined as 

Aff = D' lt Dl z \M^ z ' 1 ' + M£rrr), Mi 

Ml = D zt Dl zt (M2 zzz + M ZZlZ ), Mf 

Ml = D lt {D'p t M^ 11 + D zzt M z J Zri ), Mi 

M 7 S = D zt {D^ t M2 z ' lZ + D zzt M zzzz ), Ml 

Here. 

Df = -D il £v(k)Di l , (18) 

with (fc) being the transverse part of the renormalizcd photon, Z-boson and jZ self-energies. The longitudinal 
parts of the boson self-energy make contributions that are proportional to the ratios m 2 /t and m 2 /u so they 
are very small and are not considered here. 



= -D^ u Df\M22 Zl + M Z d 777 ), 

= -D Zu D} zt (M2 ZZZ + M z d ZyZ ), 

= -D-< u {Dj lt M22 1 ' < + D zzt M z 2 Zl ), 

= -D Zu (Dj ft M2^ Z + D zzt M zzzz ). (17) 



B. Vertices 

In order to calculate the electron vertex corrections (2nd and 3rd diagrams in Fig. [2]), we use the form factors 
SFy A in manner of paper [24| . Replacing the coupling constants v-* , a? with these form factors (for example, 

v f( z ) — >. 5Fy^ e , — > SF A ^ e ) in the corresponding terms of the Born functions Af eViOC i, we obtain a vertex 
part of the cross sections in the form 



Ver 



TTCt 

s 



1=1 



(19) 



where 



and 



Mi 



M. 



M2f r ")), 



Mo = 



M 



D^(D^{M^^ + M22 F1 ~<) - ^"(Af^ 777 + A/ 77 ^ 7 )) 
]yrt( D zt( M r>2hz + M iJF~'Z) _ D z \M F l z ^ z 

D zt (D^(M[ v Z '' z ^ + M z 22 pZl ) - D"i u (M^ lZl + M z 22 pZ ~<)), 
D zt {D zt {M F-zzz + m zzfZz ) - D Zu (M^ zzz + M ZZFZz )), 



\ v % ' 3 = 5F v vJ + SFia j , X F A h3 = SFy^ + 6F%v> . 
For the virtual photon exchange (i = 7 in Eq. (J5J)) , vertices look like 



SF 



v = 



a 
-in 



A 1 + ((v z ) +(a z Y ')A 2 (m z ) + — —A 3 (m w ) 

4s w 



(20) 
(21) 

(22) 



SFl 



a 

47T 



2v z a z A 2 (m z ) + -^-A 3 (m w ) 



and for the virtual Z-boson exchange (i = Z in Eq.([5])) we have the form 



(23) 



SF6 



a 
4tt 



v z A 1 + v z ((v z ) +3{a z ) )A 2 (m z ) + 



A 2 (m w ) - "jj^A 3 (m w ) 



is 



(24) 



SFi 



a 

47T 



a z A 1 + a z {3(v z Y + (a z ) )A 2 (m z ) + 



A 2 {m w ) - ^-A 3 (m w ) 



(25) 
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Functions Ai, the contribution of triangle diagrams with additional photon, A2, the triangle diagrams with 
additional massive boson - Z or W, and A3, the triangle diagrams with 3-boson vertex - WW"f or WW Z . 
We denote the vertices with an adititonal photon as light vertices (LV) and the vertices with an adititonal 
massive boson as heavy vertices (HV). The LV terms are proportional to the function Ai, and the HV terms 
are proportional to the combinations of functions A2 and A3 as it is evident from Eqs. (|22 [) -(|25 p . 

The contributions of separate HV, 7Z-, and ZZ- self energies are strongly depend on the details of renor- 
malization scheme and for proper account of EWC they should be taken in numerical analysis as one gauge 



independent set. For that purpose we used the on-shell renormalization scheme from Ref. [25|. It is of crusial 
importance to raise a question of dependence of EWC on details of different renormalization schemes and this 
will be addressed in our following paper. 

Next, we represent the vertex contribution as the sum of infrared (IR) divergent and finite parts using the 
identical transformation 



a 



Ver _ ^Ver _ cr Vcr^2 ^ s )) + ^er^ ^ = ^Ver + ^Ver^a ^ ^ ^ 



where A is the photon mass which regularizes the infrared divergence. Then the IR-part of vertex cross section 
will look like 



= ^Al R ^M i ° + (m U ), (27) 



2.s 



1=1 



where 



Af = Ai - Ai(A 2 s) = -2 log ± log -4, (28) 



e is the base of the natural logarithm, and 





= D yt (D jt M^' y ' r — 




Ml 


= D^(D zt M^ z - 


-D Zu M2^ z ), 


M° 


= D zt {D lt M z ^ z ~ 1 - 


- D~< u M Z £ Zl ), 


Ml 


= D zt {D zt M zzzz 


- D Zu M z f zz 



(29) 

It is worth noticing that since the Born cross section can be presented as 

,2 4 



,0 



-J2M? + (t^u) (30) 



2s 



it is easy to make sure that in the IR-part of vertex cross section the Born structure is factorized with t- and 
it-terms separated. 



C. Boxes 



In this section all formulae are written using the low energy approximation. Our numerical estimates show the 
accuracy better than 0.2% in the whole low energy region < s/s < 50 GeV. The accuracy of the approximate 
approach improves with the decrease of energy. On the contrary, the calculation of boxes using FeynArts 
and FormCalc [2l| in the region ^/s < 1 GeV suffers from a problem of numerical instability due to Landau 
singularities. In any case, for the 11 GeV relevant for the JLab experiments, the consistency of calculations for 
boxes in both approaches is obvious and discrepancy has an order of ~ 0.1%. 

Using the identical transformation, we divide the box cross section into two parts: 

_Box /'^Box _Box / \ 2 . \\ , _Box/\2 . \ _Box , „Box /oi\ 

a = (a —a (A — S-s)) + er (A — > s) = er IR + er F . (31) 
The IR-finite box cross section looks like: 

- F Box - ~ E + f E B h) + (* - «). ( 32 ) 

i=l [ij)=lk=7,Z 
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where L u = log(— s/u), the double index (ij) runs like (ij) = {1, 2, 3, 4} = {77, jZ, ZZ, WW}, and expressions 
B k . take a form 

Sf 77) = ^AI*^) + (D kt + D ku )\l k Sf 17) , 
B k hz) = D kt \ z _H\ lZ) + (D kt + D ku )X z + k Sl z) , 
B k zz) = D kt \ Bk 8\ zz) + (D kt + D ku )\f k S 2 (zz) , 

B (ww) = D *A_ S\ww) + ' + D ka )^+ k d(ww)- (33) 

The combinations of the coupling constants are given in ((5]). Let us recall the coupling constants for heavy 
boxes: 

v B = (v z f + (a z )\ a B = 2v z a z , v c = a c = l/(As 2 w ). (34) 
At s, \t\, \ u\ <C m z the corrections 6^. 2 ^ have a form: 

<5( 77 ) = L 2 s (s 2 + u 2 )/(2t) - L s u - (L 2 + n 2 )u 2 /t, 



<5( 77 ) = L 2 s s 2 /t + L x s-(Ll+n 2 )(s 2 + u 2 )/(2t), 
8\ lZ) = 8u 2 (4L (Z - I lZ ), 6 2 hz) = 8s 2 (1^-41^), 



A 1 

°(ZZ) 



where 



with 



3u 2 




3s 2 


2m 2 z ' 


°(ZZ) 


2to| 


2u 2 


X 1 


s 2 


2 ' 

m w 


°(WW) 


~ 2m 2 ' 


L s 


= log- 


- t , L x = 



Vw) = ZT2 - ' (Wff) = ' (35) 



1 1 

logp (36) 



1 Z" 1 Z" 1 1 . xz\f—u — \f]3 
Lyz = , / / ds-^log 7 =- — -= , (37) 

2V- - " Jo Jo VP xz\J—u + VP 

/3 = -ux 2 z 2 + 4(1 - z)(tz(x - 1) + m|). 

-f 7 z = i^zU-f-s- (38) 
Finally, the IR-parts of the box cross section are given by: 

v3 2 



77— box 

S 



- ^ log log ±^M<> + (39) 



v 3 



1=1 



4 



7 Z-box = a_ ^_ 8 ^ + 

2=3 

and, summing all of the IR-terms of V^-contribution, we have 

1=1 

o 4 

= io g _log_-£jW? + (t+*ti) = log — log — 5-<7°. 41 
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IV. BREMSSTRAHLUNG 

In order to get the IR-fmite result, we have to consider the diagrams with the photon emission (Fig[3]). The 
whole brcmsstrahlung cross section looks like 



a R 



4s n 2s 

i,3=~t,Z 



dv J2 (42) 



where v cut is the boundary of the region in Chew-Low diagram jl 31 ] , and 

1 r rl 3 k 

= - / T^(( fc i +n - fc 2 - k) 2 - m 2 )[M%] (43) 

is an integral over phase-space of an emitted photon with the 4-momcntum k. The squared matrix elements 
corresponding to diagrams in Fig|3]can be presented as 

Ml = {Ml - M»)(M£ - M**)+, (44) 

where minus sign is caused by identity of final electrons and the superscripts t(u) denote the i(u)-channel 
diagrams. 

Simplifying (|44|) . we have 

M%M% + = £ Vl j , M*M% + = £ V?, (45) 

fc=l,4 fe=5,8 

M$Mg + = M*M% + \ k2 « P2 , M»M£ + = Mf il4" + | fe2+ ^ 2 , (46) 

where 

V[ j = -MGrP ll {ki)GT T Kk2)]Ml^P l] ipi)lMP2)]D^D^, 
V? = -Sp[GrP y '(fci)7,A(fc 2 )]Sp[7 Al p y '(pi)Gr T A(p 2 )] J D ltl J D'''*, 
V; j = -Sp[GrP y '(Pi)7,A(p 2 )]Sp[ 7A1 ^(fc 1 )Gr T A(fc 2 )] J D lt ^* 1 , 
V? = -Sp[GrP y '(pi)Gr T A(p 2 )]Sp[ 7A1 ^(fc 1 ) 7 ,A(fc 2 )] J D' ;t J D^, 
V* = -S P [G^^(k 1 )G^A(p 2 ) 1 ^( Pl )^A(k2)]D^D^, 
V« = -Sp[Gr^(fci)7,A(p 2 )7 Al ^'bi)GrA(fc 2 )]^ 1 ^^, 
V^' = -Sp[7 A1 ^(fci)GrA(p 2 )Gr^'bi)7,A(fc 2 )]^^", 

V? = -Sp[^(k 1 )^A(p 2 )G^ a p^(p 1 )GrA(k 2 )]D it D^. (47) 

Wc used the radiative invariants which are zero at k — > 0: 

zi = 2kk\, z = 2kk 2 , vi = 2kp 1: v = 2kp 2 = s + u + t— Am 2 , (48) 

and the invariants for propagator structure as: 

<i = (pi-p 2 ) 2 , z 2 = {k 1 -p 2 ) 2 . (49) 

Also, wc applied short abbreviation for combinations with p{jp), density matrix of particle with 4-momcntum p, 
defined as 

p*t(p) = K - oS)p(p)(^ + a 3 75) = ^(Ai'jp - \4 l5 p) + 0(m). (50) 

We define 

A(p)=p + m, p = ^p fl , (51) 

and 

; 2fcf-fc 7 Q 2k%+-f a k 



= ™j + -2 / ^y , 

— Zl z 



G 2 =7 1 7 ! 

— Vl v 

G 3 = — 7 +7 , 

2pf-y4 2fc 2 " + fc 7 " 

G 4 = 7 +7 . 

— Vi z 
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n 





(3) 



FIG. 3: Bremsstrahlung diagrams for the M0ller process in the t-channel. The four u-channel diagrams are obtained from 
those in Fig. 3 by means of the interchange ft2 f>p2- 



The bremsstrahlung cross section can be broken down into two parts (soft and hard) as 



'IR 



(52) 



by separating the integration domain according to ko < oj or fco > uj, where fco is the photon energy (in the 
reference frame co-moving with the center of mass of primary electrons) and w is a parameter corresponding to 
the maximum soft photon energy. The easiest way to implement such partition is to multiply the integrand in 
(|42|) by 9(u — ko) and neglect photon momentum k — > where possible, which would give us the soft photon 
cross section. 



A. Soft Photons and IR-divergence Cancellation 



First, we follow the methods of paper |26( to get a well-known result (see also [11|, [12|) for the soft photon 
cross section: 

,2 



B a , 2lu tu 

0IR = - ( 4 log -T- log — - l0j 



.2 



7r A em^s em 

Next, we summ the IR-terms of V- and ^-contributions (formulae (|4Tj) and (|53p) 



(53) 



a Auj 2 tu 
-(2 log log 



°IR + °IR — V" *~t> 2 

7r s em^s 



loe 



1 



t*"" , 2 u \ n 

y +los i> - 



(54) 



and get a result free from IR-divergence which logarithmically depends on oj and contains log (s/m 2 )-terms. 
Adding the contribution corresponding to Ai to a c , we get an expression with the first power of collinear 
logarithms: 



^Ver , C 

a + a 



;(Ai(A 2 -> S )-log 2 ^) + 



7T 

a —t —ts 
- log — log — 

7T S 771 



with non-physical dependencies cancelled analytically. 



a /i 2 — * , 2 s x 
- log 2 — - log 2 — + 
7r m z 



(55) 



B. Hard Photons. Leading Logarithms Approach 



Now we will calculate the hard bremsstrahlung cross section retaining in the result the leading collinear 
logarithms. This approach allows estimating the EWC very rapidly yet provide a rather accurate result. First 
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let us consider exact collinear kinematics, with all relevant information given in Table 1 of Appendix A. As one 
can see from this table, z + v ~ (1 — ri)s for all peaks. Using the radiative invariants fco = (v + z)/(2y / s), we 
can transform the region of integration over v into 

1 - 2-^= < 77 < 1 - 2^=, (56) 
V s V s 

where f2 is the maximal energy of emitted hard photon. 

Now, integrating over photon phase-space and taking into consideration the results from Table 1 to simplify, 
we get the hard-part of bremsstrahlung cross section as 



Jh 4s J 1 - r/ 

l-20/\/s 



Mv), (57) 



where 



h(rj) = h CVtt (v) + h od) t(v) + ( fc 2 ^ P2), (58) 



and 

s — v( (1 - rj) 2 



h eY M= S — !L \ (2im-2 + ^ 5Lf a ) V M««( W «)i5 <tl i^ tl |* I + 

+ (2l m - 2 + (1 ~ ^ * a ) V M^(s, u/ V )D^D^ \ z 
J! MH ij (ris,u)(D Ul D jt + D lt D jtl )\ Zl + 
+(!* - L) E MH ij (s,u/ V )(D^D jt + D lt D^)\ z + 



-2( (1 2 ^ V +VQ E M c r( S ,u)^^T 

i-j=l-Z j 



K*M = V —^- E M of'( ('7('? 2 -ry+l)f a + ?? 2 ; m -2ry 2 )^ tl J D^| Zl 4 

«,3=7j 



E M of f C^-^ + l^ + r? 2 



+7j 2 (Z m - J.)^** 1 !^*!* 1 + (lx + l m + (1 L ~ 2)D^D^\ Z 



+r, 2 (l a - l s )D lt D^ + (l u - -)D lt D^ u \ v 



+((t ? 2 - V + l)l u - 2 V )D lt D^ r + (!x ~ la)D lt D^ \ z . (59) 



The logarithms above look like 



(s — v) 2 ~ (s + 1) 2 —t 

la = log 5 , l a = log 5 , l m = log — , 

l s = log 4, lu = log = log ^ (60) 



x 10 



0.12 



0.1 



0.08 



0.06 
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FIG. 4: Born asymmetry versus 9 at Ei^ = 3.03 GeV (solid line), 11 GeV (dashed line). 



with t = v + to 2 . An operation E\ x denotes the calculation of an i?-expression in the a;-peak. 
Using the standard designation of IR-divergence extracting operation (see a pioneer paper [2] 



we can present the hard part as 



R _ R,U . R,u _ a 



3 r A m ^,3 



4s 



dy , . a" 
(1-77)+ 4s 



drj 



Mi). 



(61) 



(62) 



Obviously, the fist term in (|62l) is independent of u (at sufficiently small uj it depends on Q only), and the 
second term can be easily calculated as: 



1-2 



R.u 



■Ml) 



dij ia CI 



tu 



= — log -log— ^-cr°. 
4s ./ 1 — 77 7r w cm s 



(63) 



1-2- 



Finally, we obtain the desired result where the sum of all contributions to lowest order radiative corrections 
is independent of uj: 



a c + o- 



H 



(64) 



V. NUMERICAL ANALYSIS 



Fig.|4]shows the Born asymmetry A® versus 9 at the energy of previous JLab experiment at Ey^ = 3.03 GeV 
(solid line) and the energy of planned one with E\ & \, = 11 GeV (dashed line). Here, we used a = 1/137.035999, 
my/ = 80.398 GeV, and mz = 91.1876 GeV according to (28|. It is clear that at low energies the asymmetry 
A° is proportional to s = 2mE uh giving A? (11 GcV)/,4°(3.03 GeV) w 4 for any 6. 

For detailed numerical calculation of EWC, we take electron, muon, and r-lepton masses as m e = 
0.510998910 MeV, m M = 0.105658367 GeV, to t = 1.77684 GeV and quark masses for loop contributions as 



m u = 0.06983 GeV, m c = 1.2 GeV, m t = 174 GeV, m d = 0.06984 GeV, m & 
The light quark masses provide Aq;^ (J (to 2 j)=0. 02757 [2!|, where 



u.d.s.c.b 



0.15 GeV, and m b = 4.6 GeV. 



(65) 
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Qf is the electric charge of fcrmion / in proton's charge units q, (q = y/Awa). We believe that the use of the 
light quarks masses as parameters regulated by the hadron vacuum polarization is a better choice in this case. 
Finally, for the mass of the Higg boson, we take run = 115 GeV. Although this mass is still to be determined 
experimentally, the dependence of EWC from mu is rather weak. 
Let us define the relative corrections to the Born cross section as 



S c = {<t c -<i a )/a Q , C = BSE,Ver,Box, 



and to the Born asymmetry as 



S C A 



{A°-A°)/A°. 



For the convenience, we define " weafe" for relative corrections as all BSE contributions (including 77-SE which 
is not weak by nature, but we need it here to account for all IR- finite contribution to asymmetry), HV, ZZ- 
and WW-boxes. This way, weak = BSE+UV+ZZ+WW. 

We present all weak and total (total = weak+QED) relative corrections to the unpolarized cross section at 
Eiab = 11 GeV and different angles 6 with different cuts on soft ("S") and hard ("S+H") photon emission energy 
given by various 71 (71 = Ej/y/s) in Table 1 of Appendix B. Empty cells correspond to forbidden kincmatical 
region of this process. If E 1 = u>, we take into consideration only soft photons and if E 1 — il, both soft and 



hard photons are taken into account. In principle, for _E 7 = fi, we should treat co/y/s as a small parameter but 



as it was shown in Eq. (|64p . the total result docs not depend on it in any case. The 77-SE contribution is small 
but still dominates the relative weak correction to the unpolarized cross section. We can also see a rather small 
difference between contributions "S" and "S+H". Additionally, we can conclude that unpolarized cross section 
significantly drops with the decrease of 71. 

At this point, it is essential to compare the corrected parity-violating asymmetry, which is sensitive to input 
parameters and calculation scheme, with well-known existing results. In Fig. 5 we can see the relative weak 
(solid line) and QED including soft (dashed line) corrections to the Born asymmetry A® versus yfs at 9 = 90°. 
In the region of high energies {yfs > 50 GeV) we can see an excellent agreement with the result of Dcnner and 
Pozzorini [TTJ] if we used their Standard Model parameters (see Table below for S^ eak at different yfs) . 



yfs, GeV 


result of Rcf. [11] 


our result 


100 
500 
2000 


-0.2787 
-0.3407 
-0.9056 


-0.2790 
-0.3406 
-0.9066 



Furthermore, the relative QED correction (see Fig. 8 in Ref. [ll| and dashed line in Fig. 5 here) are also in 
a good qualitative and numerical agreement. In this case, we apply the same cut on the soft photon emission 
energy as in Rcf. [llj] (7! = 0.05). 

At low energy point corresponding to E-158 experiment and using our parameters set we find that S^ eak = ~ 
-54%. If we translate our input parameters to the scheme by Marciano and Sirling according to [ll[, we obtain 
a good agreement with the result of (loj . 

All weak and QED relative corrections to asymmetry A® at 8 = 90° and Ei a b = 11 GeV are presented in 
Tables 2 and 3 in Appendix B. The 77-SE gives ~ —0.007 to the relative weak correction. The 7Z-SE gives 
considerable contribution and the ZZ-SE gives a contribution which is non-zero but not as large as from 7Z-SE. 
It is unphysical to separate the BSE and HV due to the fact that only their sum gives a gauge invariant set. 
We used the separation to tune up our code by comparing the values obtained with different renormalization 
schemes [24| and [25j]. We show their contribution as a total weak correction in Table 2. The Wiy-boxes give a 
rather significant input: S^ w = 0.0238, but on the contrary the result coming from ZZ-boxes is rather small 
and equal to <5 ZZ = —0.0013. Finally, we can see from Table 3 that the total relative QED corrections are 
significant and strongly depend on the parameter 71. 



VI. CONCLUSION 



M0ller scattering is a very clean process that can provide indirect access to new physics at multi-TeV scales. 
The new ultra-precise measurement of the weak mixing angle via 11 GeV M0ller scattering to start soon at JLab 
will require the higher order effects to be taken into account with the highest precision as well. In this work, we 
calculate the electroweak radiative corrections to asymmetry of polarized M0ller scattering at energies relevant 
to the future experiments at JLab. The results are presented in both numerical and analytical form. As one 
can clearly see from our numerical data, at certain kinematic conditions EWC can reduce the asymmetry up to 
~ 70%, and they depend quite strongly on the experimental cuts. In the on-shell renormalization scheme we 




-1 2 3 

10 1 10 10 10 



FIG. 5: The relative weak (solid line) and QED (dashed line) corrections to the Born asymmetry A® versus y/s at 6 = 
9CP . The filled circle corresponds to our predictions to the future 11 GeV M0ller experiment at JLab 



use, the largest contribution is coming from the weak corrections: self-energy graphs (especially 7Z-SE), HV, 
and boxes. The QED part is obviously very important as well. 

We believe that one of the most important result of this work, however, is our compact analytical and 
sufficiently accurate expressions. They can be useful for fast yet relatively precise estimations and are well 
suitable for building Monte Carlo generators. Our final results are analytically free from any non-physical 
parameters. The accuracy was controlled by comparison with the numerical data obtained by semi-automatic 
approach using FeynArts and FormCalc. These base languages were already successfully employed in the similar 
projects (see [30 ] and |3lj . for example), so we arc highly confident in their reliability. 

Since the problem of EWC for M0ller scattering is rather involved, we believe that the tuned step-by-step 
comparison between different calculation approaches is essential. In the next work, we plan to present a 
detailed comparison between several calculation approaches with different renormalization schemes. We also 
plan to address the leading two-loop electroweak corrections which are very likely to be required by the promised 
experimental precision. To maximize the precision, the full set of one-loop EWC evaluated in this paper will 
also need to be re-calculated with the latest input parameters available at the time of the completion of the 
measurements, but this will be easy to do based on the results of the present work. 
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Appendix A: Collinear kinematics 

Tabic 1. Kinematic rcationships relevant to hard photon cross section. 





zi-peak 


z-peak 


t>i-peak 


u-peak 


k 


(1 - rj)ki 




(i - v)pi 






2(1 - 7?)m 2 -> 


^(2m 2 -t) 


(1 — 7])(s — 2m 2 ) 


(l-r ? )(s + t-4m 2 + 


z 


(l-?y)(2m 2 -f) 


2^1m 2 -> 


±-Z( vs + t-2m 2 ) 


(1 - r])(s - 2m 2 ) 


V\ 


(1-?7)(s-2to 2 ) 


±-H( v ( s -2m 2 )+t) 


2(1 - T])m 2 -> 


(1 - ?7)(2to 2 - t) 


V 


(l-r])(s + t-2m 2 ) 


(l-r?)(s-2m 2 ) 


i^(2m 2 -t) 


2^m 2 -> 
n 2 


u 


rj(2m 2 -s-t) + 2m 2 


2m 2 - t - r](s - 2m 2 ) 


2m 2 s+ 2m ^ t 
u 


2m 2 - s - t + 

n 


k2P2 


V ■ hpi 


rj ■ kipi 


V ■ k\p\ 


rj ■ kipi 


P\P2 


r\ ■ kik 2 




1 ■ fcifc 2 


r\ ■ kik 2 


hP2 






r]-pik 2 


r\-p\ki 
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Appendix B: Numerical analysis 

Tabic 1. The unpolarized Born cross section and the relative weak and total corrections to it at £i a b=ll GeV 
at different 71 (71 = 0.005, 0.01, 0.05) and 9. 



a 


cr , mb 


weak 


Q nnr; 
b, U.UUo 


n 1 TT r\ (\C\K 

b+ri, U.UUo 


Q m 
0, U.U1 


G 1 u m 
b+ri, U.U1 


0, U.Uo 


C 1 tr n r,r 
o+ri, U.Uo 


20 


0.1277 x 10 2 


0.0087 


-0.2149 


-0.2148 


-0.1754 


-0.1758 






30 


0.2607 x 10 1 


0.0101 


-0.2417 


-0.2415 


-0.1972 


-0.1978 






40 


0.8734 


0.0111 


-0.2595 


-0.2591 


-0.2118 


-0.2124 


-0.1012 


-0.1067 


50 


0.3920 


0.0119 


-0.2721 


-0.2716 


-0.2222 


-0.2227 


-0.1063 


-0.1136 


60 


0.2176 


0.0126 


-0.2810 


-0.2805 


-0.2295 


-0.2303 


-0.1099 


-0.1183 


70 


0.1444 


0.0131 


-0.2870 


-0.2867 


-0.2344 


-0.2356 


-0.1124 


-0.1219 


80 


0.1131 


0.0135 


-0.2905 


-0.2904 


-0.2373 


-0.2389 


-0.1139 


-0.1241 


90 


0.1043 


0.0136 


-0.2916 


-0.2916 


-0.2383 


-0.2400 


-0.1144 


-0.1249 



Table 2. The Born asymmetry A® and the structure of relative weak corrections to it at £ , i a b=H GeV at 
different 0. 



e° 


A?, PP b 


77- SE 


7Z-SE 


ZZ-SE 


all BSE 


ZZ-box 


WW-box 


HV 


weak 


20 


6.63 


-0.0043 


-0.2919 


-0.0105 


-0.2863 


-0.0013 


0.0239 


-0.2946 


-0.5643 


30 


15.19 


-0.0049 


-0.2916 


-0.0105 


-0.3065 


-0.0013 


0.0238 


-0.2633 


-0.5430 


40 


27.45 


-0.0054 


-0.2914 


-0.0105 


-0.3051 


-0.0013 


0.0238 


-0.2727 


-0.5508 


50 


43.05 


-0.0058 


-0.2912 


-0.0105 


-0.3042 


-0.0013 


0.0239 


-0.2703 


-0.5489 


60 


60.69 


-0.0062 


-0.2911 


-0.0105 


-0.3045 


-0.0013 


0.0239 


-0.2714 


-0.5500 


70 


77.68 


-0.0064 


-0.2910 


-0.0105 


-0.3040 


-0.0013 


0.0238 


-0.2712 


-0.5495 


80 


90.28 


-0.0066 


-0.2909 


-0.0105 


-0.3040 


-0.0013 


0.0238 


-0.2711 


-0.5493 


90 


94.97 


-0.0067 


-0.2909 


-0.0105 


-0.3043 


-0.0013 


0.0238 


-0.2710 


-0.5493 



Tabic 3. The Born asymmetry and the QED corrections to it at £i a b=ll GeV at different 71 (71 = 
0.005, 0.01, 0.05) and 9. 



6° 


A\, ppb 


S, 0.005 


S+H, 0.005 


S, 0.01 


S+H, 0.01 


S, 0.05 


S+H, 0.05 


20 


6.63 


-0.0710 


-0.0649 


-0.0676 


-0.0566 






30 


15.19 


-0.0758 


-0.0736 


-0.0716 


-0.0686 






40 


27.45 


-0.0792 


-0.0790 


-0.0744 


-0.0744 


-0.0651 


-0.0567 


50 


43.05 


-0.0817 


-0.0826 


-0.0763 


-0.0778 


-0.0663 


-0.0660 


60 


60.69 


-0.0833 


-0.0848 


-0.0777 


-0.0797 


-0.0671 


-0.0691 


70 


77.68 


-0.0844 


-0.0859 


-0.0785 


-0.0805 


-0.0675 


-0.0701 


80 


90.28 


-0.0849 


-0.0863 


-0.0789 


-0.0806 


-0.0677 


-0.0702 


90 


94.97 


-0.0850 


-0.0863 


-0.0790 


-0.0806 


-0.0678 


-0.0702 



